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ACCESSIBILITY PERCOLATION AND FIRST-PASSAGE SITE 
PERCOLATION ON THE UNORIENTED BINARY HYPERCUBE 

ANDERS MARTINSSON 


Abstract. Inspired by biological evolution, we consider the following so- 
called accessibility percolation problem: The vertices of the unoriented n- 
dimensional binary hypercube are assigned independent U (0,1) weights, re¬ 
ferred to as fitnesses. A path is considered accessible if fitnesses are strictly 
increasing along it. We prove that the probability that the global fitness max¬ 
imum is accessible from the all zeroes vertex converges to 1 — 2 In (2 + \/E) 
as n —> oo. Moreover, we prove that if one conditions on the location of the 
fitness maximum being v , then provided v is not too close to the all zeroes 
vertex in Hamming distance, the probability that v is accessible converges to 
a function of this distance divided by n as n —> oo. This resolves a conjecture 
by Berestycki, Brunet and Shi in almost full generality. 

As a second result we show that, for any graph, accessibility percolation 
can equivalently be formulated in terms of first-passage site percolation. This 
connection is of particular importance for the study of accessibility percolation 
on trees. 


1. Introduction 


A number of recent papers [4 TO have studied a percolation problem known 
as accessibility percolation, based on ideas of Kauffman and Levin for modeling 
biological evolution [I]. In its simplest form, accessibility percolation consists of a 
graph G = (V,E), or more generally a digraph, together with a fitness function 
w : V —> R. generated according to some random distribution. This is thought 
of as representing the landscape of possible evolutionary trajectories of a species. 
The vertices in G represent the possible genotypes for an organism whose fitness 
is a measure of how successful an individual of that genotype is, and the edges the 
possible ways the genome can change subject to a single mutation. Here it makes 
sense both to consider directed and undirected edges depending on whether or not 
a certain mutation is reversible. Of primary concern is the existence or distribution 
of so-called accessible paths. 


Definition. Let G = ( V , E ) and w : V —» K. be a fitness landscape. We say that a 
path vq —> v\ —)■ ... —tvi in G is accessible if 

(1.1) w(u 0 ) < w(ui) < ... < io(vi). 

For v, w £ V we say that w is accessible from v if there exists an accessible path 
from v to w. 


For the distribution of to we will in this paper consider two variations of King¬ 
man’s House-of-Cards model |B|. Both of which have previously been considered in 
accessibility percolation. In fact, all results in [6jjl0 consider some variation of the 
House-of-Cards model, whereas and :5| also consider the so-called Rough Mount 
Fuji model. The first model we will consider here is the original formulation of the 
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House-of-Cards model, in which the w(u):s are independent and 17(0, l)-distributed 
for all v £ V. Kauffman and Levin refers to this as an uncorrelated landscape. 
For the second distribution we modify the House-of-Cards model by introducing 
an a priori global fitness maximum v £ V by changing oj(v) to one. As accessi¬ 
bility percolation only depends on the relative order of fitnesses, this can be seen 
as equivalent to conditioning the House-of-Cards model on v being the global fit¬ 
ness maximum. In particular, if v is chosen uniformly at random among V, then 
this is equivalent to the House-of-Cards model with v denoting the global fitness 
maximum. 

Our first main result considers accessibility percolation on the unoriented n- 
dimensional binary hypercube. The question of primary concern is whether or not 
there exists an accessible path from the all zeroes vertex, 0, to the fitness maximum 
v. We prove that, provided v is not too close to 0 in Hamming distance, the 
probability that such path exists converges to a non-trivial function of the Hamming 
distance between v and 0 divided by ?r, confirming a conjecture by Berestycki, 
Brunet and Shi |7j in almost full generality. 

As a second result, we show that accessibility percolation for a general graph can 
be equivalently formulated in terms of first-passage site percolation. This lets us re¬ 
formulate previous results in the literature in terms of first-passage site percolation. 
In particular, this relation has important implications for accessibility percolation 
on trees, as studied in [6j[8|jl0]. 

1.1. Notation. 

• Whenever talking about a general graph G = (V,E), we allow both undi¬ 
rected and directed edges. For vertices u, v £ V, we write u ~ v if there is 
either an undirected edge between u and v or a directed edge going from u 
to v. 

• The unoriented n-dimensional binary hypercube, denoted by Q„, is the 
graph whose vertices are the binary ?r-tuplets {0,1}” and where two vertices 
share an edge if their Hamming distance is one. The oriented n-dimensional 
binary hypercube, Q„, is the directed graph obtained by directing each edge 
in Q n towards the vertex with more ones. 

• For a vertex v in the hypercube we let v denote the number of coordinates 
of v that are one. Addition and subtraction of vertices in Q n denotes 
coordinate-wise addition/subtraction modulo two. We let 0 and 1 denote 
the all zeroes and all ones vertices respectively, and let e\ ,..., e n denote 
the standard basis. 

• Often when considering the House-of-Cards model, it is useful to condition 
on the fitness of 6. Following the convention in 6,7, for any a £ [0,1] we let 
P Q (-) and E Q [•] denote conditional probability and expectation respectively, 
given w(0) = a. 

1.2. Recent work. Let us take a moment to summarize the results for accessibility 
percolation on the binary hypercube with House-of-Cards fitnesses in §0- We 
start by consider the simplified version of the problem where we replace Q ra by Q ra . 
This is equivalent to only considering paths without backwards mutations. As any 
coordinate where v is zero will be constantly zero along any such path, it suffices 
to consider the case where v = 1. 

Let A' denote the number of oriented paths from 0 to 1 which are accessible. As 
there are n\ oriented paths from 0 to 1, and each path is accessible if and only if 
the n random fitnesses along the path are in ascending order, we see that EX = 1. 
At first glance, this may seem to imply a positive probability of accessible paths 
existing. However, a much clearer picture of what occurs is obtained by conditioning 
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on the fitness of the starting vertex. Indeed, conditioned on the fitness of 6 being 
a £ [0,1], we have 

(1.2) E a X = n(l-a) n ~ 1 . 

We see that, for large n, this expression is 1 approximately at a = ^Ji, and rapidly 
decreasing as a increases. Informally, this means that unless the fitness of the start¬ 
ing vertex is below 1,1 ”, accessible paths are highly unlikely. In fact, by considering 
( |1.2[ ) a bit more closely it follows that P (X > 1 A w(0) > < A On the other 

hand, the regime where a is smaller than ln " turns out to be more difficult to treat. 
In 5| it was shown by Hegarty and the author that the probability of accessible 
paths in this case tends to 1 as n —> oo. 

Theorem 1.1. (Hegarty, Martinsson) For any sequence {Sn}^! such that ne n —> 
oo, as n —> oo we have 

(1.3) PTr+ E »(X > 1) -> 0 

(1.4) P^-^iX > 1) -a 1. 

Furthermore, 

(1.5) P(X > 1) ~ —. 

n 

This theorem was later strengthened by Berestycki, Brunet and Shi in [6] who 
proved that, in the special case where w(0) = O(^), X has a non-trivial limit 
distribution when scaled appropriately. 

Let us now switch back to the unoriented hypercube and see how this analysis 
changes. Again, let X denote the number of accessible paths from 0 to v. Here, 
paths are not as combinatorially well-behaved as for the oriented cube, and first 
moment estimates are not as easy to come by. Nevertheless, in a recent paper by 
Berestycki, Brunet and Shi [7] it was shown that E Q X has the following asymptotic 
behavior: 


Theorem 1.2. (Berestycki, Brunet, Shi) Let a £ [0,1] be fixed, and let v = v n £ 
Q n be such that x := hin^^oo \v n \ /n exists. We have that as n —> oo 

(1.6) (E Q X) 1/n -A- sinh(l - a) x cosh(l - aY~ x . 

As a consequence, for each x there is a critical value a*(x) = l — d(x) for the fitness 
of 0 , given by the unique non-negative solution to 

(1.7) (sinhi?)^ (cosh'd) 1- ' 1 = 1, 
such that 

• For a > 1 — i9(x), P“ (X > 1) goes to 0 exponentially fast as n —>• oo. 

• For a < 1 — i9(x), E Q X diverges exponentially fast as n —> oo. 

Hence, the unconditioned probability that X > 1 is at most 1 — d(x). 


We see a similar behavior of E a X as for the oriented cube. One important differ¬ 
ence though is that unlike the oriented cube the critical value has a nontrivial limit 
as n —> oo. The function 'd(x) is plotted in Figure [I] This function is continuous 
and increasing where $(0) = 0 and i?(l) = lnp + v^) ~ 0.88. In particular, it 
follows that if the a priori global fitness maximum is 1, then the critical fitness is 
1 —ln (l + -s/2) ~ 0.12, and if chosen uniformly at random then |D| /n will be tightly 
concentrated around | and hence the critical fitness is 1 — | ln (2 + y/5) 

Berestycki et al. further gave two conjectures that (1.61 
sense that P“ (X > 1) tends to 1 as n 


0.28. 

tells the truth” in the 
oo for a < 1 — d(x). Conjecture 1 of their 
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Figure 1. The function d(x) as defined in (1.7). 


paper proposes this in the special case where •0 = 1, and Conjecture 2 in the more 
general setting of v = v n satisfying \v n \ /n x £ [0,1]. 

1.3. Results. The first result of this paper fully resolves Conjecture 1 by Beresty- 
cki, Brunet and Shi |7 , and Conjecture 2 under the additional condition that x is 
not too small. 

Theorem 1.3. Let v = v n £ Q„ be a sequence of vertices such thatx := linin^oo |t)| /n 
exists. Let X denote the number of accessible paths from 0 to v. Let i9(x) be as 
defined in Theorem [Ol Assuming x > 0.002, we have 


( 1 . 8 ) 


lim P' 

n—>oo 


(X > 1) = 


if a > 1 — x ) 
if a < 1 — d(x). 


In particular, if v = 1, then 
(1.9) P(X>1) 


-> 1 - In (l + V2 


and if v is chosen uniformly at random, then 

(i.io) - 1 


P (X > 1) -> 1 - i In (2 + V5 


The value 0.002 deserves some explanation. In the proof of Theorem |1.3[ or more 
accurately the proof of Theorem 1.6 below which is shown to be equivalent to the 
former, we see that there is a value x* « 0.00167 such that the proof goes through 
whenever x > x* and breaks down when x < x*, see Remark |4.8| It seems likely 
however that this is simply an artifact of the technique used in the proof, and that 
the statement should hold true even for smaller x. Regardless of whether or not 
this is true, we can note that the two cases of most concern, x = 1 and x = 0.5, are 
far above x*. 

We now turn to the relation between accessibility percolation and first-passage 
site percolation for a general graph. Let G = (V, E) be a graph with a distinguished 
vertex 0. Note that each edge of G may either be directed or undirected. For each 
vertex v £ G randomly assign a cost, denoted by c(v ), according to independent 
1/(0,1) random variables. For a path uq, u±, ..., ui in G we define the site passage 
time of the path by 

(1.11) C M> 

l<i<l 
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Figure 2. Example of a graph where accessible paths have a dif¬ 
ferent distribution than paths with small passage time. We can 
for instance note that there can never be exactly three accessible 
paths from 0 to v, whereas there can certainly be exactly three 
paths with reduced passage time at most 1 — a. 


and similarly define its reduced site passage time by 

( 1 . 12 ) c ( u i)- 

1 <i<l 

Note that neither the passage time nor the reduced passage time of a path include 
the cost of the first vertex. For each u,v £ G we define the site first-passage time 
from u to v, denoted by 7 v(u,v), and the reduced first passage time from u to v, 
denoted by Ty(u, v), as the minimum of the respective quantity over all paths from 
u to v. 

Theorem 1.4. Let G be a graph with two distinct vertices 0 and v, and let a £ 
[0,1]. Consider accessibility percolation on G. If fitnesses are assigned according to 
the House-of-Cards model with v as the a priori global fitness maximum, then 

(1.13) P“ (v accessible from 6) = P (7y(0, v) < 1 — a) . 

If fitnesses are assigned according to the House-of-Cards model without an a priori 
global maximum, then for any vertex v £ G 

(1.14) P a (y accessible from 6) = P (Tv (6, v) < 1 — a) . 

Moreover, in the latter case this claim can be significantly strengthened. Conditioned 
on the fitness of 0 being a, the set of vertices accessible from 0 has the same 
distribution as the set of vertices v such that 7y ( 0 , v) < 1 — a. 

Informally we can think of this theorem as saying that accessibility percolation is 
equivalent to first-passage site percolation with independent 1/(0,1) vertex passage 
times. We need to be a bit careful there though; the theorem only deals with the 
question of whether or not a certain vertex is accessible from 0 along any path, 
and it does not for instance say anything about the number of accessible paths. 
Indeed, it is not true in general that the number of accessible paths from 0 to v is 
distributed as the number of paths from 0 to v with reduced passage time at most 
1 — a. For graphs containing non-simple paths this is clear as non-simple paths 
can have arbitrarily small passage time but cannot be accessible, but it can even 
be false for directed acyclic graphs, see for instance Figure [2] On the other hand, 
the connection is more general than just treating which vertices are accessible. For 
instance, using the proof ideas in Section [2] one can show that the minimal number 
of times you need to move to a less fit vertex to get from 0 to v is distributed as 
the integer part of Tv(0, v ) + a. 

A problem with using Theorem |1.4| to relate known results from first-passage 
percolation to accessibility percolation is that the vast majority of the first-passage 
percolation literature assigns passage times to edges rather than vertices. However, 
a common property for percolation problems is that it is harder to percolate on 
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vertices than edges (2|. The following proposition shows that something similar 
holds for first-passage percolation. 

Proposition 1.5. Suppose the edges of G are assigned independent 17(0,1) weights. 
Let Te(u, v ) denote the minimum total weight of any path from u to v in G. Then, 
it is possible to couple Te(0 , •) to 7y(0,-) such that Te(0, v) < Tv(0,v) for all 
v G G. 

For the special case when G is a rooted tree one can see that this coupling 
is exact; to go from site to bond percolation we can simply consider the passage 
time of each vertex to instead be assigned to the edge leading to it. Accessibility 
percolation on trees has been considered in |:6][8j]l0 . With the exception of [6], 
these articles have considered regular rooted trees with degree n and height h, and 
where fitnesses are assigned according to the House-of-Cards model conditioned on 
the fitness of the root being zero. Of principal concern is how the number of vertices 
in generation h that are accessible from the root varies as a function of n, and in 


particular whether this number is non-zero. Using Theorem 1.4 we can see that this 
is equivalent to assigning independent 17(0,1) passage times to the edges of the tree 
and considering the number of vertices v in generation h such that Te(0, v) < 1. In 
particular, the question of whether generation h is accessible from 0 is equivalent 
to to asking if the first-passage time from the root to generation h is at most 1. It 
should be mentioned however that the usual setting in first-passage percolation on 
regular rooted trees keeps n fixed and considers the first-passage time from the root 
to generation h as h —> oo. While the author is not aware of any results from this 
field that have appropriate error bounds to be directly applicable to accessibility 
percolation, there seems to be a significant overlap of ideas between 9 10 and the 
literature on first-passage percolation on trees. See for instance 13]. 

Let us now consider the implications of Theorem |1.4| for the hypercube. Using 
this result, we can immediately translate the result from Theorem o to that, for 
the oriented hypercube, Ty (0,1) is concentrated around 1 — 7m w ith fluctuations 
of order -. More importantly, we have that the following is equivalent to Theorem 

773 


Theorem 1.6. Let G = Q n and let v = v n £ Q n be a sequence of vertices such 
that x := lim n _). 00 |D| /n exists. Assuming x > 0.002, as n oo we have 


(1.15) 


Ty(0,v) i9(x) 


in probability. 

Note here that the fact that i9(x) is an asymptotic lower bound on the reduced 
passage time is already implied by Theorem |1.2[ 

It should be mentioned that basically the same results holds true for bond perco¬ 
lation. In 11 it was shown that for the oriented hypercube, we have Te(0, 1) —)• 1 
in probability as n — > oo. In a more recent result by the author 12], it was shown 
that for the unoriented hypercube 7 e( 6, 1) —» In (l + \/2) as n — > oo. Strictly 
speaking these results assume standard exponential edge weights, but it is not too 
hard to show that the limiting distribution of Te(0, 1) only depends on the weight 
distribution as the rightlrand limit of its probability distribution function at 0, hence 
it will be the same for U( 0,1) weights. 

The remainder of the paper will be structured as follows: In Section [2] we prove 
Proposition 1.5 and Theorem 1.4 The remaining sections, Sections [ 7 ] [4]and[5j are 


dedicated to the proof of Theorem 1.6 
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2. Proof of Proposition 11.51 and Theorem 11.41 

We may, without loss of generality, assume that for any vertex v there exists a 
path from 0 to v. 

A key idea of the proofs of Proposition |1.5| and Theorem |1.4| is the following 
procedure for computing 7v(0,v). We initially consider Tv(0,v) to be unassigned 
for each v, except 0 for which it is set to 0, and we let U = {0} denote the set of 
vertices with assigned first-passage times. Until Tv{Q,v) is assigned for all v , we 
do the following operation: 

(1) Find a pair of vertices u,v that minimizes Tv(0,u) subject to {u,v) £ E, 
u £ U and v £ U. 

(2) Let Tv (6, v) := Tv (6, u ) + c( v) 

(3) Add v to U. 

To see that this assigns first-passage times correctly, suppose that we are in the step 
where 7y(0, v) is assigned. As v is not in 17, the passage time of any path from 0 to v 
must include the passage time from 0 to some vertex u' in U adjacent to some vertex 
outside 17, as well as the cost v. Hence 7y(0,u) > Tv(0,u')+c(v) > 7v(0, u)+c(v). 
As there is a path from 0 to v with passage time 7y(0,u) + c(v), this must be 
optimal. Hence, if all previous assignments are correct, 7y(0,v) will be assigned 
correctly as well. 

Proof of Proposition \ 1.5\ We can modify this algorithm to run on first-passage 
bond percolation by replacing c(v) by the weight of the edge from u to v. In 
either case, as no vertex cost or edge weight respectively is accessed more than 
once, the accessed values form a sequence of independent and 17(0,1) random vari¬ 
ables. Hence the distribution of Tv(0, v ) is unaffected. On the other hand, for bond 
percolation we get that Tv(0, v) is the edge passage time of some path from 0 to v 
(but not necessarily the shortest). ■ 

We now turn to the proof of Theorem |1.4| The coupling between first-passage 
site percolation and accessibility percolation we will consider is essentially to let 
f(v ) = {a + 7y(0,u)} be the fitness function, where {.t} = x — \x\ denotes the 
fractional part of x. We will however modify this slightly by putting f{y) = 1 
whenever a + 7y(0,u) = 1. It is clear that the probability of such v other than 0 
existing is 0, so the only way this will change the distribution of / is that /(0) = 1 
if a = 1. 

It is not too hard to see that, for any vertex v except 0, f(v) is U (0, l)-distributed. 
The following lemma shows that the f(v): s are also independent, hence showing that 
/ is distributed according to the House-of-Cards model without an a priori global 
fitness maximum, conditioned on /(0) = a. 

Lemma 2.1. f(v) are independent U( 0,1) random variables for v £ V \ {0}. 

Proof. Suppose that we generate vertex costs in the following way: Run the pro¬ 
cedure above, but with the modification that whenever the algorithm tries to ac¬ 
cess c(v), first generate a 17(0,1) random variable f(v) and assign c[v) the value 
{/(u)-a-7V(6,u)}. 

It is clear that the c(v):s are independent and 17(0, l)-distributed. The lemma 
follows by noting that, in the latter case, we have f(v ) = f(v ) almost surely for all 
w£V\{0}. ■ 

Proof of Theorem \l.f\ We begin by considering the case with no a priori global 
fitness maximum. In this case, we can consider / : V -I R to be the fitness 
function. For simplicity let us assume that no vertex cost is exactly 0. 
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Assume 7 r(0, v) < 1 — a, and let 6 = Vo, iq ,..., Vi = v be the path with shortest 
passage time. Then, as 0 < a + 7v(0, Vi) < 1 for 1 < i < l it follows that 

i 

(2.1) f(v i ) = a + ^2c(v j ) 

i=i 

for 0 < i < l. Hence vq , Vi, ... , Vi is accessible. Conversely, suppose Tv{0, v) > 1 — a 
and let 0 = vq , Vi, ... ,v i = v be any path between 0 and v. Let i be the lowest index 
such that Tv(0,Vi) > 1 — a. Then /(iq) < a + Ty(0,Vi) — 1 < a + Ty(0,t; i _ 1 ) + 
c(vi) — 1 < a + Tv (6, Vi- 1 ) = f{vi- 1 ). Hence the path is not accessible. 

Now for the case where v £ V \ {0} is the a priori global fitness maximum. We 
here keep the same coupling as before between f(v) and c(v ) for v € V, except that 
we fix f(v) = 1. Let U be the set of vertices v £ V such that (v,v) £ E. Then 
v being accessible from 0 is almost surely equivalent to some vertex in U being 
accessible from 0. Note that this last statement does not depend on the value 
of f(v). It follows that v is accessible from 0 is almost surely equivalent to that 
min^giy Tv(0, v) < 1 — a. The theorem follows by noting that min„g[/ 7y(0,v) = 
6,15). ■ 


3. The Clustering Translation Process 


Before proceeding, we will slightly modify T(((6,f)) by replacing the 17(0,1) ver¬ 
tex costs by independent standard exponential such. Note that the standard expo¬ 
nential distribution stochastically dominates 17(0,1), and h ence this modification 
will only increase Ty{0,v). As the lower bound in Theorem 1.6 follows from The¬ 
orem |1.2[ it suffices to show that, with this modification, asymptotically almost 
surely Ty(0,v) < i9(x) + o(l). To do this, we will mimic the argument in 12j for 
first-passage bond percolation on Q n . 

Let us take a moment to describe some of the underlying machinery for first- 
passage bond percolation on Q n . We assume independent standard exponential 
edge weights. In 11 , Durrett introduced the following process, which he called 
the the branching translation process, BTP: At time 0 we place one particle at 0 
in Q n . The system then evolves by each existing particle independently generating 
offspring at each vertex adjacent to its position at rate 1. One can show that for 
each vertex v £ Q n , the time at which the first particle at v is born is stochastically 
dominated by Te(0,v). This follows from the fact that the BTP dominates the 


so-called Richardson’s model. The strategy in 12 is basically to show that, with a 


certain coupling, there is a probability bounded away from zero of these quantities 
being equal. 

In order to translate this approach to first-passage site percolation, we need to 
find a corresponding process to the BTP for this case. We claim that the following 
is such a process: We initially have a finite number of particles, each located at a 
vertex in Q n . For each particle, we assign an independent Poisson clock with unit 
rate. When a particle’s clock goes off, it simultaneously generates one new offspring 
at each vertex adjacent to its position. The new particles are then assigned new 
Poisson clocks and the process continues. We will refer to this process as the 
clustering translation process, CTP. 

We see that in both the BTP and CTP each particle generates offspring at each 
neighboring vertex at rate 1. A big difference however is that in the BTP this 
is done independently for each neighboring vertex, whereas in the CTP a particle 
generates offspring all neighboring vertices simultaneously. Another difference is 
that the initial state of the CTP is not fixed. 
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The most important initial state of the CTP will be one particle at each neighbor 
of 0. We will refer to a CTP initialized in this way as a standard CTP. Particles 
born due to the same Poisson clock tick will be referred to as identical n-tuplets. To 
simplify terminology we will also consider the initial n particles in a standard CTP 
as identical n-tuplets. Below we will use the terms ancestor and descendant of a 
particle to denote the natural partial order of particles generated by the CTP. For 
convenience, we say that a particle is both an ancestor and a descendant of itself. 
The terms parent and child are defined in the natural way. The ancestral line of 
a particle x is the ordered set of ancestors of x, and we say that the ancestral line 
of x follows the path 0 = Vq,v\,V 2 , ■ ■ ■ ,vi if the location of the ancestors of x in 
chronological order is given by v±, i> 2 , ■ ■ ■ Vi. Note that this path always starts at 0 
even though the first ancestor is located at a neighbor of 0. We say that a particle 
x originates from a particle y at a time t if y is the last particle in the ancestral 
line of x that exists at time t. 

We can immediately note some properties of this process. Firstly, it is Markovian. 
Secondly, let A be a set of vertices in Q n , and let MA(v,t) denote the expected 
number of particles at vertex v at time t > 0 in the CTP initialized by placing one 
particle at each vertex in A. Then it is easy to see that M A (v : t) must solve the 
initial value problem 

(3.1) = £ M A (w,t) for t > 0 

w~v 

(3.2) M a (v,0) = 1 a { v ). 

In particular, if A = {0}, then the unique solution to this problem is 

(3.3) m(v,t) := (sinht)^ (coshf)" - ^ , 
and it follows by linearity that for any A, we have 

(3.4) M A (v,t ) = m(v — w,t). 

w^A 

Recall that addition/subtraction of vertices in Q n are interpreted as coordinate- 
wise addition/subtraction modulo 2. It should be remarked that the exact same 
analysis holds for the BTP. 

We now show that the standard CTP indeed has the desired relation to first- 
passage site percolation. To this end, we partition the particles in this process into 
two sets, the set of alive particles and the set of ghosts. Each initial particle is alive. 
Whenever a new particle is born, it is alive if its location does not already contain 
an alive particle and its parent is alive, and is a ghost otherwise. Note that at most 
one particle at each vertex can be alive. Furthermore, it is easy to show that each 
vertex will almost surely eventually contain an alive particle. 

Proposition 3.1. Consider first-passage site percolation on Q n with exponentially 
distributed costs with unit mean. It is possible to couple this process to the standard 
CTP such that for each vertex v except 0, Ty(0,v) denotes the birth time of the 
alive particle at v. 

Proof. For each vertex v, we let 7~'(v) denote the first time t > 0 when v contains 
an alive particle, and we let c(v) denote the time from the birth of this particle to 
the first arrival of its clock. Then c(v) for v £ Q„ are independent exponentially 
distributed random variables with unit mean. 

From the definitions of the CTP and alive particles, it follows that for any vertex 
v that is not a neighbor of 0, the alive particle at v is born at the first arrival time 
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of an alive particle at an adjacent vertex. Hence, for any v that is not a neighbor 
of 0, we have 

(3.5) T'(v) = min (t'(w) + c(w) \ , 

and trivially T’(v) = 0 when v is a neighbor of 0. It is easy to see that this uniquely 
defines T'(v), and that for each vertex v except 0, T'(v) denotes the reduced first- 
passage time from 0 to v with respect to the vertex costs given by c(v). ■ 

Given this proposition, we are able to proceed analogously to Sections 2 and 3 in 
12 . In applying this coupling between the CTP and first-passage site percolation 
we will consider a stronger and more tractable property than aliveness. For any 
particle x in the CTP, we let c(x) denote the number of pairs of particles y and z 
such that 

• y and 2 occupy the same vertex 

• y is an ancestor of x 

• y was born after z. 

We furthermore let a(x ) denote the number of such pairs where z is either an 
ancestor of x or an identical n-tuple of an ancestor of x , and define b(x) = c(x) — 
a(x). We call a particle x uncontested if c(x) = 0. 

It can be noted that a(x) is defined differently for the BTP. This is because 
the strategy is loosely speaking to let a(x) denote the number of pairs (y, z) that 
deterministically must exist given x. For the CTP we have additional such pairs, 
namely those corresponding to identical ro-tuplets of ancestors of x. 

Lemma 3.2. If a particle is uncontested, then it is alive. 

Proof. If a particle a; is a ghost, then it must have an earliest ancestor (possibly 
itself) which is a ghost, y. As y is a ghost but the parent of y is alive, it follows 
that the location of y must have already been occupied by some (alive) particle 2 . 
The pair (y, z) is then counted in c{x). I 

It is not hard to see that a(x) only depends on the path followed by the ancestral 
line of a;. If we know this path, then we know the locations and order of births of all 
ancestors and identical n-tuplets of ancestors of x. Let a be a path represented as a 
vertex sequence. We say that er is vertex-minimal if there is no proper subsequence 
which is a path with the same end points. 

Lemma 3.3. Let x be a particle in the CTP. If the ancestral line of x is vertex- 
minimal, then a(x ) = 0. The converse is true unless x is located at 0. 

Proof. Denote the path followed by the ancestral line of a: by i) 0 , iq ,... ,i’i and the 
ancestors of x by X\, x-i,... ,xi = x. We have that a(x) > 0 if and only if there exist 
1 < * < j < l such that Xj occupies the same vertex as either aq or an identical 
n-tuplet of Xi, that is, Vi -1 and Vj are adjacent. Hence, if a(x) > 0 the path is not 
vertex-minimal. Conversely, if a(x) = 0 it follows that the only pairs of adjacent 
vertices are consecutive in the path. It is straight-forward to show that, unless the 
path starts and stops at the same vertex, this implies vertex-minimality. ■ 

What follows are two technical lemmas, corresponding to Lemmas 3.2 and 3.3 in 
12 . Before presenting these, we need to specify how to formally describe the CTP. 
Firstly, by a (potential) particle we mean a word {tq, Zi, V 2 , 22 ,... zi-\,vi) 

where tq ,... ,Vi denote vertices and 21 ,... 2 ;_ 1 positive real numbers. This is in¬ 
terpreted as the particle whose ancestors are located at v\,V 2 ,... ,vi and born at 
times 0, 21,21 + 22 and so on. The CTP is described by a random set X of poten¬ 
tial particles, denoting the set of particles that will ever be born in the CTP. We 
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will use ® to denote concatenation of words. We remark that this representation 
means that the functions c(x ) and b(x) are not functions only of x, and should more 
correctly be denoted by c(X, a;) and 6(X, x). On the other hand, a(x) is really a 
function of x as it only depends on the location of the ancestors of x. 

Lemma 3.4. Let a = {0 = vo, Vi, ..., u/_i, u;} be a path. For 0 < i < l — 1 let X, 
denote independent CTP:s where Xj is the CTP obtained by initially placing one 
particle at each neighbor of Vi . Let f be a function that maps pairs ( X , ir) to the 
non-negative real numbers where X is a realization of a CTP, and x is a particle 
in X. Similarly, let V tT (X) denote the set of particles in X whose ancestral lines 
follow cr. Then for a standard CTP, X, we have 



where 


(3.7) X‘ 1,, " Zi 1 — X 0 U ({viZi} ® Xr) U • • • U ({ V 1 Z 1 V 2 Z 2 ■ ■ ■ vi-izi-i} ® X;_iJ 

and x 21 ’--’* 1 - 1 = {V 1 Z 1 V 2 Z 2 ■ ■ ■ vi}. 


For compactness, we will only sketch a proof. The reader unconvinced by this is 


referred to the proof of Lemma 3.2 in 12 


Proof sketch. Let us first consider the case when /( X,x) only depends on x. In 
that case, we have 


(3.8) 



/ (ar* 1 ’-’*'- 1 ) dzi... dzi- 1 . 


This is because the original particle at v± gives birth to particles at V 2 at rate one 
whereupon, after its birth, each child at V 2 of this original particle gives birth to 
particles at V 3 at rate one, and so on. When / also depends on the realization of 
the CTP, the idea is that we substitute /(X, x) in the left-hand side of this sum by 
E [f(X.,x)\x £ X]. Now, formally this conditioning does not really make sense, but 
its meaning is intuitively clear; it denotes the average value of f(X,x) where the 
average is taken over all X that include x. We have 


E 22 f( X ’ x )= E E E[/(X,*)|*eX] 
x6V„(X) xev„(x) 


= ... E [/(X, x 




)k 


2l,...,2l_l 


£ X] dz\ ... dzi- 1 . 


Jo Jo 

Now, x exists in X if and only if certain Poisson clocks have arrivals at 
certain times. By the independent increment property, conditioning on these ar¬ 
rivals does not affect the Poisson clocks at any other times. Hence, the conditional 
distribution of X given the existence of x Zl >- -’ Zl - 1 is the same as that of a standard 
CTP, except with added arrivals, corresponding to the births of the ancestors of 
. This is precisely the distribution of X 21 ’"'^- 1 . ■ 


Lemma 3.5. Let X be a CTP, and let <fi be an indicator function on the set of 
potential particles in X. If cj>{x) = 0 for all original particles in the CTP, then 

(3.9) P(£kk = o)>exp(-E5>0r) 

VxGX / V xGX 

Proof. Let us refer to the set of original particles as generation one, their children 
as generation two and so on. Let T denote the set of birth times for particles in 
generation two in X, and let T' C T be the subset obtained by including t £ T 
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if there exists a particle cc £ X such that <j>(x) = 1 and x is an descendant of a 
particle in generation two born at time t. It is clear that |T'| < and 

that |T'j =0 if and only if </H a: ) = 0. 

By definition of the CTP, it is clear that T is a Poisson point process. Further¬ 
more, as the event that t £ T is included in T' only depends on descendants of 
particles in generation two born at time t, this occurs independently for each t £ T. 
Hence, by the random selection property, T' is also a Poisson point process. This 
implies that 


(3.10) P j 53 00) = 0 J = P (T' = 0) = exp (-E |T'|) > exp j -E 53 00) ) , 
VxGX / V i6X / 


as desired. 


Theorem 3.6. Consider a standard CTP. For any vertex v and any t > 0, let 
B(y,t) = E b(x) where the sum goes over all particles at v at time t in the 
CTP, and let S(v,t ) denote the expected number of particles x at v at time t such 
that a(x) = 0. The probability that there is an uncontested particle at v at time t is 
at least S(v, t) exp ^ . 

Proof. Let P(v,t) denote the probability that v contains an uncontested particle at 
time t. As at most one particle at each vertex can be uncontested, this is the same 
thing as the expected number of uncontested particles at v at time t. For each path 
a from 0 to v , let P a (v,t ), B a (v,t) and S a (v,t) denote the contribution to P(v,t), 
B{v,t) and S(v,t) respectively from particles whose ancestral line follows a. 

The idea now is to bound P a (v, t) in terms of B a (v, t) and S a (v, t) for each path 
a from 0 to v. Recall that a{x) is constant over all particles x whose ancestral line 
follows a fixed a. We will denote this constant by a(a). 

Let a be a path from 0 to v such that a(a) = 0. Applying Lemma 
that 


3.4 


we see 


(3.11) 


/*00 fOO 

Ida (v, t) — E ^ / -^SiH- \-zi — 1 dZ\ ... dzi— 1, 

ie^(x' 


where T(x) denotes the time of birth of x. Similarly, for any x £ V^X) we have 
Ba{v,t) = E 53 1 T(x)<tb(X,x) 

(3.12) zeVRx) 

v 7 n oo /*oo 

= / ■■■ l, 1+ ... + , 1 _ 1 < t E6(X* 1 --* | - l ,®* 1 '->*«- 1 ) d Zl ... dzi-! 

Jo Jo 

and 

(3.13) 


Pa(v,t) = E 


xev„(x) 


lT(x)<tlfc(X,x)=0 


; P (b (X- 




) = 0) dz!... dzi _ 1 . 


/ * ’ * / ^- 21-1 - \-Zl — 1 <£ ^ 

Jo Jo 

As a(x zi ’" ’ zi ~ 1 ) = 0, no two ancestors of x zi, "'’ zi ~ 1 occupy the same vertex. It 
follows that any pair of particles y and z which is counted in b ) 

is uniquely determined by z. Fixing a and zi ,..., zj_i, this means that we can 
define <j>(x) as an indicator function such that 


b (X 




) = 


E 

lEX" 1 .*1-1 


(f{x). 


(3.14) 
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More precisely, (f>{x) is the indicator for x occupying the same vertex as an ancestor 
of and being born before it. By the definition of b(x zi ’"' ,zi ~ 1 ), we have 

that 4 >{x) is zero for any ancestor or identical n-tuplet of an ancestor of 
It follows by Lemma [375] that 

(3.15) P(b(X Zl ’-’ Zi -\x Zl -’ Z1 - 1 ) = 0) > exp (—E6 (x zi> —>* 1 - 1 , X Z1 ’ "’ ZI - X )). 


By convexity of the exponential function we have e r > (1 + rg — r)e r ° for any 
r, ro £ R. Hence 

(3.16) P<r{v,t) > (1 + r 0 )e~ r °S a (v,t) - e~ r °B a (v, t), 


for any path a from 6 to v such that a(a) = 0. For any a that satisfies a(a) ^ 0 
it is clear that P a {v,t) = S a {v,t) = 0 and B a (v,t) > 0, hence (3.16) holds in this 
case as well. Summing over all paths a from 0 to v and optimizing over yields 
P(v,t) > S(v, t)e~ s ( v ’V , as desired. ■ 


We will apply Theorem 3.6 as follows: Let {v n }^ =1 be a sequence of vertices 
such that, for each n, v n £ Q n and x = limn-^ |D„| jn exists and is non-zero. We 
may, without loss of generality, assume that v n is never equal to 0. For each n, we 
let 0 n denote the unique non-negative solution to 

1 


(3.17) 


m(v n , 0 n ) = 


n 


Note that the expected number of particles at v n at time 0 n in a standard CTP 
on Q n is 0(1), and that 0 n —» 0(x) as n —>■ oo. By Theorem 3.6 we have that the 


probability that there is a uncontested particle at v n at t ime 0 n in the CTP on 


is at least S(v n , 0 n ) exp ( — ) ■ Hence by Lemma 

follows that 


3.2 


and Proposition 


3.1 


it 


(3.18) 


P (Ty(0,D„) < i9 n ) > S(v n , i9 n ) exp ( - 


B(v ni ‘On ) 

S{v n ,0n) 


This means that if we can show that S(v n ,O n ) = 0(1) and B(v n , 0 n ) = 0(1), then 
we know that 7^(0, v n ) < 0 n with probability bounded away from Oasn-> oo. 

Sectio n [4] w ill be dedicated to estimating S(v n , 0 n ) and B(v n , 0 n ). The proof of 
Theorem 11.6| is then completed in Section^lby showing that if Ty(0, v n ) < 0 n with 
probability bounded away from 0, then a slightly larger upper bound on Ty{ 0,v n ) 
must hold asymptotically almost surely. 


4. Calculus 

4.1. Estimating S. We will prove that S(v n ,0 n ) = 0(1) in two steps. Firstly, we 
show that most particles at v n at time 0 n have ancestral lines which are close to 
vertex-minimal. Using this, we then give a combinatorial argument that shows that 
a positive proportion of these particles must have vertex-minimal ancestral lines. 

Let us formalize the notion of paths being close to vertex-minimal. Let v, w £ Q n 
be fixed distinct vertices and let er = {v = Vq,v\, .. . ,vi = w} be a path from v to 
w. Throughout this section, we will always think of a path as a finite sequence of 
vertices. In particular, by the length of a path we mean the number of vertices in 
the path. For any 0 < i < j < l we say that the subsequence Vi, 1 ^+ 1 , •.., Vj is a 
detour of a if removing these elements from a results in a valid path. Clearly, for 
t) / w a path is vertex-minimal if and only if it has no detours. Inspired by this, 
we say that a path is almost vertex-minimal if all detours have length at most 2. 
Note that as Q n is bipartite, any detour must have even length. Hence, a path is 
almost vertex-minimal if it only has the shortest possible detours. 

An important property of almost vertex-minimal paths is that any such path 
from v to w can be constructed by taking a vertex-minimal path with the same 
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end-points and extending it as follows: Between each two adjacent elements in the 
sequence either do nothing or insert a detour of length 2. 

Lemma 4.1. Let s,t > 0 and v £ Q n . Then 



(4.1) 


Proof. Fix s. Observe that equality holds when t = 0 and that both expressions 


solves (3.1) 


Proposition 4.2. Let {{) n }^Li he a sequence of vertices, i> n £ Q n , such that 
a = limn^oo \v n \ /n exists and is positive. Then, as n —> oo, the expected number 
of particles in the standard CTP on Q n which are at v n at time d n , hut that do not 
have almost vertex-minimal ancestral lines tends to 0. 

Proof. Let X n denote the number of triples of particles x, y, 2 in the CTP on Q n 
such that 

• ai is at v n at time d n 

• y and z are located at adjacent vertices 

• z is an ancestor of y which is an ancestor of x. 

• y and z are neither one nor three generations apart. 

We note that if the ancestral line of a given particle x at v n at time can be 
constructed using some detour of length d > 2, then it is clear that x would have 
a pair of ancestors at adjacent vertices which are d + 1 generations apart. This 
means that any such x is counted at least once in X n . Hence, it suffices to show 
that EX n = o(l). 

For each triple x, y, z as above there are uniquely defined particles c, the particle 
after z in the ancestral line of x, and p, the parent of y. Note that the require¬ 
ment that y is neither the child, nor the grand-grandchild of z implies that p is a 
descendant of c, but not a child of c. 

Let T = {0 = to < t\ < • • • < tk = $n} denote the end-points of a partition 
of [0, i? n ) into left-closed right-open subintervals, and let X„ t x denote the number 
of triples as above where c and y are the only ancestors of x born during their 
respective time intervals. Pick o, b integers between 0 and k — 1. Consider the 
number of triples counted in X n> T where c is born during [t a ,t a + i) and y is born 
during [tb,tb+i). Note that this is trivially 0 whenever b < a. 

Let us count the expected number of corresponding triples for a < b. As z and 
y are located at adjacent vertices, for each such triple we may denote the locations 
of z, y, c and p by v, v + e*, v + ej and v + e* — e*, respectively for some v £ Q n 
and 1 < i,j,k < n. A particle is a potential z if it is born before time t a , hence 
there are on average X)/Li m { v ~ ei,t a ) potential z:s at v. For each z, a particle is 
a potential c if it is a child of z born during [t a ,t a + 1 ). Hence for each potential z 
at v, there are on average t a +i — t a potential c:s at v + ej. For each potential c, a 
particle is a potential p if it originates from c at time t a+ \ and is born before tb, 
but is not a child of c. Hence for each potential c at v + ej there are on average 
m(e.i — ek — ej,tb — t a+ 1 ) potential p: s at v + e* — ek if v + ej and v + ej — e*, are not 
adjacent, and m{e.i — ek — ej,tb — t a + i) — (tb — i a +i) if they are. Lastly, for each 
potential p , a particle is a potential y if it is a child of p born during [tb, tb+i), and 
for each potential y a particle is a potential x if it is located at v n , originates from 
y at time tb+i, and is born before time d n . Hence for each potential p at v + — e*, 

the expected number of potential y: s at v + e* is tb+\ — tb, and for each potential y 
at v + ei, the expected number of x:s is m(v n — v — e*, — tb+i)- Combining all 
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of these, we see that 
(4.2) 

EX„ iT = m ( u _ e; ’ ta){t a +i ~ to) ( m(ei - e k - e jt t b - t a+1 ) 

a<b vGQn i,j,k,l ' 

1 |ei+ej +efc|=1 (j'b ^a+l) j (^6+1 ^b^TYliv n ^ ^fo+l)* 


where the sums over i, j, /c and / all go from 1 to n. Letting Ti,T 2 ,... be a sequence 
of increasingly finer partitions of [0, $„] such that the length of the longest interval 
in T k tends to 0 as k —»• oo, it follows by monotone convergence that we have 
EA'„ = limfe->oo EX n t. Combining this with equation (4.2), and recognizing the 
right-hand side as a Riemann sum, we get 


EX„ = 


(4.3) 


Lemma 


4.1 


/ ^2 ^ m(v - e h a) ( mfe -e k - e j: b - a) 

0 ueQ„ i,j,k,l '■ 

- l|e i+ej+efc | = i (b - a)^m(w n - v - ei,i9 n - b) dadb. 

implies that we may replace the factor Y^veQ m(v—ei, a)m (v n — v — e*, fl r 


in the integrand of equation (4.3) by m(v n + ei + e;, d n — b + a). Hence, by the 
substitution t = b — a, the right-hand side of (4.3) simplifies to 


(4.4) 


(dn ~t) ^2 TO (®« + e i + e h d n - t)- 


■ (m(ei + ej + e k ,t) - t\ ei+ej+ek \ =1 t) dt. 
Using the fact that sinhf < cosht for all t G R, we have 


J2 m (v n + ei + e u d n - t ) (m(ei + ej + e k ,t) - l\ e . +e . +ek \ =1 t) 

i,j,k,l 

< n (sinh(i?„ - t)) 1 "" 1-2 (cosh(tf n - t)) n ~ |l ’" l+2 ^ (m(ej + ej + e k ,t) - t\ e . +e . +ek \ =1 t) . 

i,j,k 

It is straight-forward (but messy) to show that JA . k ( m(ei + e-j + e k , t) — l| e4 ^_ ei+efc | =1 t) 
(cosh£) n O ( n 3 t 3 ). As cosh (id n — t ) coshf < cosh$ n it follows that 

(4.5) EX„ < / n (sinh($„ — t) cosh t)^ v d~ 2 (coshi? ra ) n_ ^"' +2 O (n 3 t 3 ) dt. 

Jo 

Recall that by the definition of d n we have 

(4.6) (sinhi? n )^ n ^ (cosh$ n ) n_ ' u "' = —. 

n 

Define the function f(t) = lnsinh(i? n — t) + In coshf. Note that /'(£) = — coth(i? n — 
t) + tanhf, and f"{t) = — csch 2 ($ n — t) + sech 2 t. As 0 < sechf < 1 and cschf > 1 
for all 0 < t < d n it follows that / is concave, and thus for any 0 < t < d n we have 
f(t) < /(0) — t coth$ n < /(0) — t. Hence 

(4.7) (sinh(i?„ — t) coshf)' 1 ’"' -2 < (sinhi?„)^ n ' _2 e~^ Vn ^~ 2 ^. 

Plugging this into equation ( |4.5[ ), we get 

r# n 

(4.8) EA„ < / e -(KI- 2 )‘o(n 3 t 3 ) dt. 

Jo 

As \v n \ ~ x • n this implies that EA„ = O (^), as desired. ■ 


b ) 
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Proposition 4.3. For any pair of sequences {v n }ff =1 and as above, we 

have S(v n ,d n ) = 0(1). 


Proof. Let T„ and T n denote the sets of vertex- minim al and almost vertex-minimal 
paths from 0 to v n respectively. Using Lemma 3.4 with f(X, x ) as the indicator 
function of x being born at time d n and having ancestral line in T n and T„ respec¬ 
tively, we can write the expected number particles at v n at time d n in the CTP 
whose ancestral lines are almost vertex-minimal as 


(4.9) 


E 

<ref„ 


,qM-2 

On 

(M-2)! 


and the expected number that are vertex-minimal as 


(4.10) 


> 1-2 




As the total expected number of particles at v n at time d n in the CTP is 0(1), 
Proposition 4.2 implies that the sum in (4.9) is also 0(1). 

according to which 


The idea now is to group the terms of the sum in ( |4.9| 
vertex-minimal path a it is an extension of, that is we write 


(4.11) 


E 

ere f„ 


(M-2)! 


^ E E 

o-er„ 

<50 cr 


,qUI-2 

vn 

(M-2)!‘ 


Here a 3 a denotes that a is an extension of a. Note that the inequality comes 
from the fact that er may be an extension of more than one vertex-minimal path. 

Let us fix a vertex-minimal path a G T n consisting of l vertices. It is straight¬ 
forward to show that the number of possible detours of length 2 that can be inserted 
between each adjacent pair of elements in a is 3 (n — 1). Hence, there are at most 
3 fc (n — l) fc ( l k 1 ) ways to extend a to an almost vertex-minimal path of length l + 2k. 
This means that 


v- n o-i\ e 2+2fc 

Ep^E 3 ("- 1 ) 


d-er. 

oOc 


k =0 


/-I 




k ){l- 2 + 2k)\ 
l - 1\ dl k 


k =0 


< 


M, 2 A 3d 2 n (n-l) 

(l- l) 2 
3 d 2 n (n- 1) 


( 1 - 2)1 

(l-2)\ 


k ) (l- l) 2 ^ 

/-I 


exp 


/-I 


As any path from 6 to v n must have length at least |t) n | + 1, we conclude that 


(4.12) 


E 

rrer„ 


d 


M-2 


(M-2)! 


> exp —3d 


,n—l 


E 

ctG f„ 


d 


kl-2 


(M-2)! 


= ©(!)• 
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4.2. Estimating B. 


Proposition 4.4. For any v £ Q n and any u > 0 we have 
(4.13) 



ei,t) m (A — ej,t) m(v — A, u — t) dt 


U 


+ / (u — t) m (A — efc — ej,t) m( A, t)m(i) — A — ej,u — t) dt 


o 


AeQ™ i,j,k 



+ (u — t) y to(A — ei — ej,t) m {A — efc, t)m(v — A — ej,u — t) dt, 
AeQ„ 

where the sums over i, j, k and l go from 1 to n. 

Proof. We observe that B(v,u ) is bounded by the expected number of triplets of 
particles x, y, 2 in the CTP such that 

• x is at v at time u 

• y is an ancestor of x 

• y and 2 occupy the same vertex 

• z was born before y. 

Note the similarity to the quantity X n in Proposition |4.2[ For the sake of compact¬ 
ness, we will be less rigorous here, and refer to the proof of that proposition to see 
how to formalize this argument. 

Let us start by considering the number of such triples x, y and z where z has 
no ancestors in common with x and y, that is, for some i j we have that x and 
y originate from the original particle at whereas z originates from the original 
particle at e^. Denote the common location of y and z by v, and pick k such that 
the parent of y is located at v — e Note that as z is strictly older than y, y cannot 
be an original particle and hence has a parent. The lineage of x,y,z is illustrated 
in Graph 1 of Figure [3} 

Let us count the expected number of such triples corresponding to a fixed v and 
where y is born during the time interval [f, t + dt). The potential z:s corresponding 
to a fixed j are simply the descendants of the original particle at ej that are at v 
at time t. Hence the expected number of such particles is m(v — ej,t). Similarly, 
for a fixed i the expected number of potential y: s is given by m(v — eu — &i,t) dt , 
and for each potential y the expected number of potential x:s is m(v,u — t). As 
the potential z: s are born independently of the pairs of potential x:s and y:s, we 
see that the expected number of triples x , y, z that do not have common ancestors, 
corresponding to a fixed vertex v and a fixed time interval [t, t + dt) is given by 


n n n 



(4.14) 


i= 1 j =1 k—1 
j¥=i 


The total expected number of triples x , y, z without common ancestors is hence 
given by summing this expression over all vertices v £ Q n and integrating over t 
from 0 to u. This is clearly bounded from above by the first term in the right-hand 
side of equation (4.13). 

We now consider the cases where the three particles x, y, z have common ances¬ 
tors. Denote the last common ancestor of the particles by l and its location by v. 
As x and z have common ancestors but neither is a descendant of the other, there 
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Graph 3. 


Graph 4. 


Figure 3. Illustration of the possible ways x, y and 0 can be 
related. The left-most arrows describe the ancestors of z and the 
right-most the ancestors of x and y. Graph 1 shows the case when 
z has no ancestor in common with x and y. Here v is the common 
location of y and z , and v — ek is the location of the parent of y. 
For Graphs 2-4, v denotes the location of the last common ancestor 
of x and z, and w the common location of y and z. Graph 2 shows 
the case where the ancestral lines of x and z split by the birth of 
a new ancestor of x, Graph 3 the case where this occurs by a new 
ancestor of z and Graph 4 the case where the first unique ancestors 
of x and of z are born simultaneously as part of the same group of 
identical n-tuplets. 


must be a time s when the ancestral lines of x and z split. There are three possible 
ways in which this can occur, as illustrated by Graphs 2-4 in Figure [3j either a 
new ancestor of x is born, a new ancestor of z is born, or new ancestors of x and 
z are identical n-tuplets and therefore born at the same time. Observe that, in all 
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three cases, y must be born strictly after this time. We let w denote the common 
location of y and z. 

We now count the expected number of such triples corresponding to fixed vertices 
v and w, where the ancestral lines split during the time interval [s, s + dt) and such 
that y is born during [s + t, s + t + dt). The potential l :s are the particles in the 
CTP at v at time s, hence the expected number of potential l :s is rn { v ~ e i^ s )■ 
For each potential l, the probability that it gives birth during [s, s + ds) is ds. Now, 
for each possibility for the ancestral lines of x and z to split, conditioned on the 
process at time s + ds, the pairs of potential a;:s and y :s originate from a different 
particle than the potential z:s. Hence these are born independently. By following 
the ancestral lines as illustrated in Graphs 2-4 in a similar manner as above, we see 
that the expected number of triples with common ancestors corresponding to fixed 
v and w, fixed time intervals, and corresponding to each case for how the ancestral 
line splits are given by 

n n n 

(4.15) y^ y^ y^ m{v—ei, s) m{w—ek—v—ej,t) m(w—v, t) m(v—w, u—s—t) ds dt 

i—1 j —1 k—1 

n n n 

(4.16) y^ y^ y^ m{v—ei, s) m(w—ek—v, t) m(w—v—ej,t) m(v—w, u—s—t) ds dt 
i= 1 j —1 k—1 


(4.17) 

n n n n 

EEEE m(v—ei, s) m(w—ei — v — ej,t ) m(w—v—ek,t) m(v—w, u— s—t) ds dt 

i—1 j=l k—1 1=1 

k^O 


respectively. The total expected number of triples x , y, z with common ancestors is 
hence given by summing these three expressions over all pairs of vertices v, w G Q n 
and integrating over all s and t such that s, t > 0 and s + t < u. 

It only remains to simplify these expressions. We observe that summing (4.15), 

removes all dependence on s. 


(4.16) and (4.17) over all v,w € 
particular the sum of (4.15) over all v,w £ 
by summing over A = w — v and applying Lemma |4.1| we have 


Consider in 
By substituting summing over w 


(4 - 18) _ 

m(v — e-i, s) m (A — e& — &j,t ) m(A, t) m(v — A — v,u — s — t) ds dt 

L7,AGQn 

= y; to(A — efc — ej, t) m( A, t)m(v — A — ej,u — t) ds dt. 

AeQ„ i,j,k 


Integrating this expression over all s,t > 0 such that s + t < u, we see that the 
expected number of triples of particles x,y,z as above corresponding to the case 
illustrated in Graph 2 in Figure [3] is given by 

r u 

(4.19) / (u-t)J2J2 m (A — efc — ej,t) m{ A, t)m(v — A — e^, u — t) dt. 

AeQ„ i,j,k 


Proceeding analogously for (4.161 and ( 4. 17| ) we see that the expected number of 
triples corresponding to Graphs 3 and 4 in Figure [3] are given respectively by 


(4.20) 


f u 

/ (u — t) y^ m (A — efc, t) m (A — ej,t)m(v — A — g, u — t) dt 

JO A /-m : .■ 


A£Qn i,j,k 
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and 


(4.21) f (u — t) y'' y'' m(A — ei — ej, t) m(A — et, t)m(v — A — e,, u — t) dt. 

Jo A ^ , 


AeQ n i,j,k,l 


The expressions in (4.19)-(4.211 are clearly bounded from above by terms 2-4 re¬ 
spectively in the right-hand side of equation (4.131. M 


Consider the sum m (A, a ) 2 rn(v — IS., b). For any v £ Q„ we let v l denote 

the i:th coordinate of v. Define the function mi : {0,1} xK —> R. by mi(0, t) =cosht 
and mi(l,t) = sinhi. Using the fact that m(v,t) = n^=i we see that 

y~~] m(A, a) 2 m(v — A, b) 
agQ„ 

n 

= ^ mi (A®, a) 2 m\{v l + A 1 , b) 
agQ„»=1 

n 1 

= ^ mi (6, a) 2 mi (t) 1 + 6 , b) 

i—1 <5=0 

= (cosh(a) 2 sinh(&) + sinh(a) 2 cosh (b)^ k (sinh(a) 2 sinh(fr) + cosh(a) 2 cosh(6))" k 
= e nb Q cosh 2a- ie“ 26 ^ ^ cosh2a + ^e" 2 ^ 

where k = |D|. Let 

(4.22) 

G x (a,b) = a;In (cosh(a) 2 sinh(6) +sinh(a) 2 cosh(&)) 

+ (1 — x) In (sinh(a) 2 sinh(fe) + cosh(a) 2 cosh(6)) 

= b + a; In [ - cosh 2a — -e _2b ) + (1 — x) In ( - cosh 2a + -e _2& 

\2 2 J 2 

Then 

(4.23) yy m(A, a) 2 m{v — A, b) = exp (nGk (a, b)j . 

AgQ„ 

Proposition 4.5. For any e > 0 there exists a constant C E > 0 only depending on 
e such that whenever u £ [e, 1] we have 

(4.24) 

B(v, u)<C e [(n 4 t 3 + n 3 t + n 2 ) (u — t) + ( n 3 t 3 + n 2 t + n)] exp (nGk (t, u — t)\ dt , 


where k = 14)1 


Proof. The idea of the proof is to use equation (4.23) to reformulate equation (4.13) 
in terms of partial derivatives of G x {a, b). Note that by the fact that m{v, t) satisfies 
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© we have 

m (A — efe — ej, a) m(A, a)m(v — A — a, b) 

AeQ n i,j,k 

+ yy m(A — efc, a) m(A — ej,a)m(v — A — e,, 6) 

AeQn ij,fc 

= ^ m"(A, a) m(A, a)m!(v — A, 6) + m'(A, a) m'(A, a)m'(v 


AeQ„ 


1 9 3 

2 9a 2 <% 


E 

AeQn 


m( A, a) 2 m(v — 


A, 6) 


1 9 3 

= 2^6 eXP 


(nG k (a, 6)^ 


A, 6) 


Similarly, using the fact that all derivatives of m(v,t) are non-negative, we have 
E] E m (^ ~ ej — Cj, a) m(A — e*,, a)m(v — A — e*, 6) 

AeQn ij.fe,; 


i a 4 

< 


6 da 3 db 


exp (riG k (a, b) 


and 


E/ E TO A — efc — e*, a) m(A — ej,a)m(v — A, 6) 

A£Q n i,j t k 

1 9 3 / \ 

- e^ 3 exp r G ^ (a,6) J ■ 

Let c denote the minimum of \ cosh 2a — | e~ 2b over all a, b > 0 such that 
e < a + b < 1. It is clear that c > 0. This means that for any a, b in this range and 
any 0 < x < 1 we have 
(4.25) 


l G 


sinh 2a 


1 cosh 2a — he 2b 


i„- 2 b r 


sinh 2 a 


cosh 2a + ie -26 


< c 1 sinh 2a. 


Hence, for sufficiently large C > 0 we have \-^G x (a,b)\ <C a whenever 0 < x < 1 
and a,b > 0 such that e < a + b < 1. Moreover, as G x (a,b) is smooth wherever it 
is defined, we know that for C sufficiently large all partial derivatives of order up 
to 4 of G x (a, b ) are bounded in absolute value by C when the pair (a, b) is in this 
domain. 

By explicitly writing out the partial derivatives of exp (^nG k(a,b)^j above and 
combining this with Proposition 4.4 we see that (4.24) holds for sufficiently large 
C, as desired. ■ 


For a given sequence v = v n as above, we define 
fn{t) = G k (t, $ n - t) 


(4.26) 


= i? n — t 4— In ( - cosh 21-e 

n 

n — k 


-20„+2t 


n 


In ( i cosh 2 t+\e~ 2 ^ +2t 
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and 


(4.27) 


f(t) = G x (t,&(x) - t) 


= d(x) — t + xln 
+ (1 - x) In ( 


^cosh2f- I e “ 2l?(x ) +2 ^ 

-cosh2i+-e- 2 ^+2A 
2 2 


Note that / depends on x. From the definition of G x (a , b) we see that /„(0) = — ^ 
and fni'&n) = -2^, and that /(0) = /(i?(x)) = 0, see ( fTT] ). 

Suppose that / n (f) is “asymptotically U-shaped” in the sense that exists a con¬ 
stant A > 0 such that for sufficiently large n we have 


(4.28) f n (t) < max (/ n (0) - At, /„($„) - A(i9„ - t)) 

for any 0 <t <‘d n . If this holds, then by Proposition |4.5| we have 

/»00 

(4.29) B(v n ,i} n ) < / O (n 3 f 3 -f n 2 t + n) e~ Xnt dt, 

Jo 

which would imply that B(v ni fin) = 0(1) as desired. It remains to show for which 
sequences of vertices v = v n , /„ is asymptotically U-shaped. We start by giving a 
simple sufficient condition for x. 


Proposition 4.6. Suppose that x > 1 — ~ 0.054. Then /„(f) is asymptoti¬ 

cally U-shaped. 


Proof. By some straight-forward but tedious calculations we see that 


(4.30) 


/n(*) = -- 


4(1 - 2e -2 ”' 1 ) 


-k 4(l + 2e~ 2 ^) 


n (cosh2f — e - 2 i?„+ 2 t^ 2 n (cosh2t + e -2, U+ 2t ) 2 

Now, if we assume that > In y/2, then the first term in the right-hand side is 
non-negative, and so we have that /"(f) is at least, say, for all 0 < t < It 
follows that if z?(x) > In y/2, then /„(t) is asymptotically U-shaped. The proposition 
follows by the easily verified fact that i9 (l — = l n \/2- H 


It is clear from the proof of Proposition 


4.6 


that the limit 1 — lrl j 2 ^ 2 ' 1 is not 


optimal, and can be lowered by considering /''(f) more closely. It turns out however 
that there is a limit for x at which the convexity of f n breaks down, and more 
importantly for sufficiently small x the asymptotic U-shape of f n breaks down. In 
the remaining part of this section, we will investigate when this occurs. 

By some more straight-forward but tedious calculations we see that 


fn(t) ■ (cosh2t — e 


-2i9„+2 t 


) (cosh 2 1 - 


,~2p n +2t 


(4.31) 


= - - e 


-40 r 


41 


-At 


e —-e — 


- e 

n — 2 k 

:-< 

n 


) 


— 


This expression has the same sign as /^(f). We see that depending on the sign of 
| — e -4 ”" it is either increasing or concave, hence f n changes sign at most twice. 
Furthermore, if f n changes sign twice it goes from negative to positive to negative. 
In the same way, since 


f'(t) ■ (cosh 2 1 - e - M (“)+ 2t ) (cosh 2t + e~ 2 ^ +2t ^ 


= I L _ e -4>3(x) 


(4.32) 


the same must be true for /(f). 


4t _ f -At 

4 


- - + 2(1 - 2x)e~ 2 ^ x) . 
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Figure 4. Plot of equation (4.34) divided by a; as a function of 
$(x). We see that as d tends to 0, this converges to its limit of 
—4. The curve intersects the d-axis at d(x) « 0.0898, that is at 
x « 0.00167. 


Combining this observation with the fact that /"(t) is bounded it follows that 
a necessary and sufficient condition for f n being asymptotically U-shaped is that 
lim^oo /;(0) = /'(0) < 0 and lim^oo /,)($„) = /'(i?(x)) > 0. In fact, the former 
condition is implied by the latter as then /^(f) changes sign at most once, but 

m = mx))=o. 

As 1 = (coslid) (tanhd)^ we have 

, jno , In cosh i? d 2 „( i? 4 

(4.33) x(d) = - 


— In tanh d — 2 In d 


O 


(In d)‘ 


Hence, we have an explicit expression for (4.32) as a function of d(x). Plugging 
t = d(x) into the right-hand side of this expression we get 


(4.34) 


je" - + 2(1 - 2*)e- M - 


Note that this has the same sign as f(d(x)). By Taylor expanding this expression 
in x and d we see that the dominating term for small d is —4x. Hence /„ is 
not asymptotically U-shaped for sufficiently small x. To get a picture of what 


happens when x increases, we divide (4.34) by x and plot as a function of d, see 


Figure |4j It is clear that there is a critical value x* slightly less than 0.0017 such 
f n is asymptotically U-shaped if and only if x > x*. This proves the following 
proposition: 


Proposition 4.7. Let be a sequence of vertices, v n G 

limn^oo \v n \ jn exists and is strictly greater than x*. Then for d r 
(3.17) we have B(v n ,fl n ) = 0(1). 


J n , such that 
as defined in 


Remark 4.8. Throughout this section we have only really been interested in de¬ 
riving a tractable upper bound for B(v n ,d n ) without discussing sharpness. Never¬ 
theless, it is not too hard to convince oneself that the bound given in Proposition 
4.5 is sharp up to, say, a polynomial factor in n. However, for x < x* we know that 
there exists an interval of positive length for t where f n (t) is positive, which would 
then imply that H(i) n ,'d n ) diverges exponentially fast in n. 


5. Completing the proof of Theorem 11.6 


Let be a sequence of vertices, v n G Q„ for each n, such that x = 

lim^oo \v n \ /n exists and is at least 0.002 and let be as in ( |3.17 ). Applying 

the estimates of S(v n , i9 n ) and B{v n , i9 n ) from Propositions |4.3|and|4.7| to Theorem 
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3.6 it follows by Proposition 3.1 and Lemma 3.2 that there exists a constant cq > 0 
such that 

(5.1) liminfP (Ty(0,v n ) < i9 n ) > c 0 . 

Since d n -A i9(x) as n — > oo, this means in particular that for any e > 0 we have 

(5.2) liminfP (Ty(0,v n ) < i9(x) + e) > c 0 . 

n—> oo 

Note that we can assume that cq is independent of the choice of sequence. 

Proposition 5.1. Let {t)n}^Li be a sequence as above, and let x = \im.\v\ /n. 
Then, for any £ > 0 we have 

(5.3) P (7^(6, fi„) <0(ar)+e)-H 


Proof. Let £ > 0 be arbitrary. Condition on the vertex passage times of all neighbors 
of 6 and v n . Assuming \v n \ > 3, it is easy to see that the number of coordinate 
places 1 < i < n with the property that the i:th coordinate of v n is 1, and the cost 
of both d and v n — are at most er/3, is distributed as Bin (|i)„| , (1 — e -£ / 3 ) 2 ). 
Hence as n —>■ oo it is clear that, with probability 1 — o(l), there are at least two 
such coordiantes. Pick a pair i ^ j. 

Depending on the choice of i and j, we define Qq as the induced subgraph of Q n 
with vertex set {v £ Q n : Vi = 1 ,Vj = 0}. We similarly define Q\ as the induced 
subgraph of Q„ with vertex set {v £ Q„ : Vi = 0 ,Vj = 1}. Note that Q 0 and Q\ 
are vertex disjoint subgraphs of Q n , both isomorphic to Q n _ 2 . 

In light of Qo and Q i, we have two natural upper bounds for Ty(0, v n ), namely 
c(ei) + c(v n — ej) plus the smallest reduced vertex passage time for any path from 
Ci to v n — ej in Q 0 , and c(ej) +c(v n — e^) plus the smallest reduced vertex passage 
time for any path from Bj to v n —et in Q\. As the only vertices of Qo and Q\ which 
are neighbors of 0 or v n are e,, ej, v n — and v n — ej, the reduced first-passage 
times in Qo and Q i are independent of each other and each is distributed as the 
reduced first-passage time between two vertices at distance |D n | — 2 in Q„_ 2 - By 
applying (5.2) to the first-passage percolation problems in Qo and Q i, we conclude 
that for any £ > 0 and for any sequence where v n £ Q„ for each n > 1 

such that x = linin^oo |D ra | /n exists and is at least 0.002, we have 

(5.4) 


liminfP (Tv(0,v n ) < d(x) + e) > 1 — (1 — c 0 )^ 


Note that this is the same expression as (5.2), except that the right-hand side here 
is strictly larger. Hence, by iteratively applying this argument, we see that we can 
replace the right-hand side in (5.4) by Ck = 1 — (1 — Co) 2 for any non-negative 
integer k. The Proposition follows by letting k —» oo. M 
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